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Abstract 

Contact of a fluid with a solid or an elastic wall is investigated. The wall 
exerts molecular forces on the fluid which is locally strongly nonhomogeneous. 
The problem is approached with a fluid energy of the second gradient form 
and a wall surface energy depending on the value of the fluid density at the 
contact. From the virtual work principle are obtained limit conditions taking 
into account the fluid density, its normal derivative to the wall and the curvature 
of the surface. 

1. Introduction 

The phenomenon of surface wetting is a subject of many experiments [1]. 
Such experiments have been used to determine many important properties of 
the wetting behavior for liquid on low energy surface [2]. In fact the wetting 
transition of fluids in contact with solid surfaces is an important held of re- 
search both for mechanics and physical chemistry. In the recent paper [3], the 
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first author using statistical methods proposed an explicit form for the energy 
of interaction between solid surfaces and liquids. This energy yields a bridge 
connecting statistical mechanics and continuum mechanics. To obtain boundary 
conditions between fluid and solid it is also necessary to know the behaviour of 
the fluid as well as the solid. 

We propose a mechanical model similar to that used in the mean-field theory 
of capillarity that leads to the second gradient theory of continuous media in 
fluid mechanics [4]. The theory is conceptually more straightforward than the 
Laplace one to build a model of capillarity [5,6]. That theory takes into account 
systems in which fluid interfaces are present [7]. The internal capillarity is one 
of the simplest cases since we are able to calculate the superficial tension in the 
case of thin interfaces as well as in thick ones [8]. It is possible to obtain the 
nucleation of drops and bubbles [9]. 

It seems that the approximation of the mean-field theory is too simple to be 
quantitatively accurate. However, it does provide a qualitative understanding. 
Moreover, the point of view, that the fluid in interfacial region may be treated 
as a bulk phase with a local free energy density and an additional contribu- 
tion arising from the nonuniformity which may be approximated by a gradient 
expansion truncated at the second order, is most likely to be successful and 
perhaps even quantitatively accurate near the critical point [10]. 

In this paper we connect both the interaction of a solid surface and a fluid 
phase by means of the virtual work principle. The distribution of fluid energy in 
the volume and the surface density energy on the solid surface yield boundary 
conditions. The conditions arc different from those obtained for a classical fluid 
within the theory of gas dynamics. We obtain an embedding effect for the 
density of the fluid; moreover, the conditions take into account the curvature of 
the surface. The result is extended to the case of an elastic wall. 

A discussion is obtained depending on the value of the density of the fluid at 
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the surface. 



Let us use asterisk "*" to denote conjugate (or transpose) mappings or cov- 
ectors (line vectors). For any vectors a, b we shall use the notation a*b for 
their scalar product (the line vector is multiplied by the column vector) and 
ab* or a (g) b for their tensor product (the column vector is multiplied by the 
line vector). The product of a mapping A by a vector a is denoted by A a. 
Notation b* A means covector c* defined by the rule c* = (^4* b )*. The 
divergence of a linear transformation A is the covector divA such that, for 
any constant vector a, 

div(A) a = div (A a ). 

If /(x) is a scalar function of the vector x associated with the Euler variables 
df 

in the physical space, ^ is the linear form associated with the gradient of / 
and Ik? = ^fx^' Consec L uentl y' (|~)* = § rad / 



2. Continuous mechanical model of capillary layers 



We consider a fluid in contact with a solid. The fluid occupies the domain D 
and its boundary X which is common with the solid wall. Physical experiments 
prove that the fluid is nonhomogeneous in the neighbourhood of £ [10]. It is 
also possible to consider the fluid as a continuous medium by taking into account 
a capillary layer existing in the vicinity of £ and a form of its stress tensor 
[11]. One way to present the behaviour of such a fluid is to consider the specific 
internal energy £ as a function of the density p as well as grad p. Such an 
expression is known in continuum mechanics as internal capillary energy, see 
[4,5]. It is related to molecular models of strongly non homogeneous fluids in 
the frame of the mean field theory and is equivalent to the van der Waals model 
of capillarity (see the review by Rowlinson and Widom [10]). The energy e is 
also a function of the specific entropy. In the case of isothermal media at a given 
temperature, the specific internal energy is replaced by the specific free energy. 
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In the mechanical case the entropy or the temperature are not concerned by 
the virtual variations of the medium. Consequently, for an isotropic fluid, it is 
assumed that 

e = f(p,P) 

where f3 = (grad p) 2 = gradp . gradp (the dot . denotes the scalar product). 
The fluid is submitted to external forces represented by a force potential SI per 
unit mass as a function of Eulcr variables x. 

We denote by x £ £ — ► B(x) £ R the surface density of energy of the solid 
wall. The total energy E of the fluid in D and its boundary £ is the sum 
of the three potential energies: E = Ef + E p + Es with, 

E f = J J J pefaP) dv, E p = J J J P^) dv , E s = J J^B ds 

Let us denote by 5 a variation of the position of the fluid as in [12]. The 
variation is associated with the virtual displacement 

x e D -> <5x = C{x) 

We have the following results presented in Appendix, 

SE f = [ [ [ (- div a) . C dv + 



D 



with 



//J- 4 "! 1 +(|An + gra,V4 + „, (1) 

a = -PI - C gradp ® gradp = -PI - C'A* ^ 

ax ax 

where C — 2pe'p and P = p 2 e' p — p div(C grad p), 

e denotes the partial derivative of e with respect to p, Q n = n*^ where n is 

dp dp 

the external unit normal to S and A = Cp — where — = — — n. 

dn dn ox 

The scalar R m is the mean curvature of S and grad tff is the tangential part 

of grad relatively to S. 

Moreover, 

5Ep = I J Id ^ = / 1 Id 9 gradr! ' C ^ (2) 
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5E s =JJsB-( ^n + grad ts s).C^ (3) 



and using the results presented in Appendix, 

2B 

One assumes that the volume mass in the fluid has a limit value p s at the wall 
E. One assumes also that B is a function of p s only. These hypotheses are 
confirmed by results presented in [3]. Then 

SB = B'(p s )Sp s = -p s B'(p s ) divC 

Let us denote G = —p s B' Ps . Consequently, 

JBds = J J^Gdw(d S = J J (G^- |^n. C-grad^G. C) ds 

(see Appendix). 

Now, H — B{p s ) — p s B' Ps (p s ) is the Legendre transformation of B with 
respect to p s . Then, 



(4) 



The d'Alembert-Lagrange principle of virtuals works is expressed in the form 
[12]: 

V x G D -» C(x), SE = (5) 

Consequently, from the fundamental lemma of variation calculus, we obtain the 
balance equation in the fluid D and the boundary conditions on the solid wall 
E. 

Equilibrium equations : 

From any arbitrary variation x G D — > C( x ) such that C — on we take 
first 

Consequently, 

-diver + p— =0 (6) 
ax 



This equation is the well known equilibrium equation [5,7,9] 



Boundary conditions : 

a) Case of a rigid (undeformed) wall. 

We consider a rigid wall. Consequently, the virtual displacements satisfy on E 
the condition n* £ = . Then, at the rigid wall 
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(G-A)^ + { ^J 1 ^ + ura,L,;.l - H ) + an ; . , do = 



Hence, we deduce the boundary conditions at the rigid wall 

For x G E, G - A = (7) 

and moreover, there exists a Lagrange multiplier xeS^ A(x) £ i? such that 
2{A-H) 



R r , 



n + gradj (A - H) + <rn = A n (8) 



b) Case of a elastic (non-rigid) solid wall. 

In such a case the equilibrium cquation(6) is unchanged. On S the condition (7) 
is also unchanged. The only different condition comes from the fact that we do 
not have anymore the slipping condition for the virtual displacement (n*£ = 0). 

Due to the possible deformation of the wall, the virtual work of mechanical 
stresses on E is 

SE e = 



with T = Q n representing the stress (loading) vector, where Q is the value 
of the Cauchy stress tensor of the wall on the boundary E. 



Relation (8) is remplaced by : 



R„ 



n + grad t (A- H)+an= T (9) 



G 



3. Analysis of the boundary conditions 



Relation (7) yields : 



Cf n+ B>=0 ,10, 



and we obtain, 

H - A = B. 
Consequently, from the definition of a, 

crn = Pn — C^f- gradp. 
an 

The tangential part of equation (8) is always verified and equation (8) yields 
the value of the Lagrange multiplier A. 

For an elastic (non-rigid) solid wall we obtain 

J5-0 and r„ = g + F- S ;£ ,11) 



where T£ g and T n are respectively the tangential and the normal components 

*2B 1 

of T. Taking into account equation (10), T n = P + — h — (B') 2 and 

R m C 

equations (11) yield the value of the stresses in the elastic (non rigid) medium. 
The only new condition comes from equation (10). 

We have the consequences: 

In [3] we propose the surface energy in the form B(p) = —-fip + —p with 
71, and 72 as two positive constants. We obtain the condition for the fluid 
density on the wall 

° tl = 71 ~ 129 (12) 

Denoting again p s the value of p on the wall S (the limit value of the volumic 

mass of the fluid on the wall £ ) , we obtain that — — is positive (or negative) 

dn 

in the vicinity of the wall if p s < pi (or p s > pi ) with pi = — which is the 

72 

bifurcation fluid density at the wall. 
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If p s < pi we have a lack of fluid density at the wall. If p s > pi we have a 
excess of fluid density at the wall. 

Conclusion 

For conservative medium, the first gradient theory corresponds to the case 
of compressibility. To take into account superficial effects acting between solids 
and fluids, we propose to use the model of fluids endowed with capillarity. The 
theory interprets the capillarity in a continuous way and contains Laplace's 
theory. The model corresponds for solids to elastic materials with couple stresses 
indicated by Toupin in [13]. 

We notice that the extension to the dynamic case is straightforward: by virtual 
work principle, equation (6) takes the form: 

/?7 - div cr + p — = 0, 
ox 

where 7 denotes the acceleration of the fluid. Equations (10), (11), (12) and 
consequences in paragraph 3 are unchanged. 
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Appendix 



First of all we recall the following fact issued from differential geometry: 

Let S be a surface in the 3-dimensional space and n its external normal. For 

any vector field (, 

2 dC 

n*rot(n x C) = div C + n*C - n* n. 

R m ax 

Then, for any scalar field A, we obtain: 



dC 2 A 

A div C = A - —Cn- (grad^A) C + n* rot (An x C) 
an Kn 



tr 



(Al) 



Let us calculate SEf : since D is a material volume, 

Ef = J J J p s dv => £E/ = y y y p Se dv 

ds 3s 

with fe = op 5/7 + ~dp 5l3 ' 

j dp cMp dp <9C . . 

From o — = — — — , we deduce 

ax ox ax ax 

, „ , „,ap ap* n (d5p dp d( \ dp* 

with 2ps'p = C. 

In the mean-field molecular theory, the quantity C is assumed constant [10], 
but it is not necessary for our calculations. One can suppose the scalar C is a 
general function of p and even (3. Then 

pe'p 5(3 = div(C grad p Sp) - div(C grad p)Sp-tr(^C gradp grad*p ^ ) 

Due to the fact that Sp = — pdiv £ (see [12]), 

pSe = div(C grad p Sp) — (p 2 e' p — p div(C grad p )jdiv ( 
— div(C gradp grad*p () + div(C gradp grad*p)C 
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pSe = div(c gradp Sp — (C grad p grad* p)( — P(j 
dP 

+ ( + div(C grad p grad* p) ( 



Then 



dv 



6Ef = J J i D ( fa + diY ^ C grad/9 grad * p) ) C 
J J div^Cp grad p div £ + C grad p grad*p £ + P<^ 

= J J J -{diva)(dv + J J {-A divC + nV C) da 
Taking into account (A.l), we deduce immediately 

SE f = J J J -(diver) C dv 

+ J J (-A^ + (^-n*+ grad*,, A + nV) () ds + J J n*vot(A n x ()ds 
But J J n*rot(A n x Q ds = J At . (n x Q d£ = J T A(t, n, Q dl 

where T is the line boundary of S and t its tangent unit vector. If n' = t x n 
we obtain the relation 



SE f =J J j (-divcr)C dv + 

J J^(-A^ + (^-n*+ gmd* tg A + n*a)()ds + ^A n'*C dt (A.2) 

In the following we assume that S has no boundary and consequently the term 
associated with L vanishes. 
Let us calculate 5E$ 

E s = J J Bds 

Then 

6ES = 1 1 SB ^( n *^- + S Tad * B ( 1 - nn *))C ds + J An'*(d£ (A3) 
We notice that grad* 5(1 — nn*) belong to the tangent plane to S. 
Let us proof equation (A. 3): If we write E$ = J J B det (n, d\x, d 2 x) , where 
dix and d 2 x are the coordinate lines of X, we may write, 

E s = [ [ B det F dot (F _1 n,diX,d2X) 
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where E is the image of E in a reference space in Lagrangian coordinates X 
and F the deformation gradient tensor — — . 

Then, 

SE S = I I 5B det F det (-F _1 n, d\~K, d 2 X) + 

J Je 

J J B<y(det F det (F~ 1 n,d 1 X,d 2 xfj 

Moreover, 

y J BS(detF det (F- 1 n,d 1 X,d 2 xfj = 

[ [ _B div^ det (n, c?ix, rf 2 , x) +_B det (^C, rfix, d 2 x)— -B det ( ^-n, dix, d 2 x) 
J Jt, " ox dx 

div(S C) - grad*B C - B n* ^ n ) ds 

ox / 



From (A.l) we obtain, 

div(B()-B (divn) n *C-n*^ n = n* rot (Bnx() 

Then, 

f [ B5(dct F det (F- 1 n,d 1 X,d 2 X)) = 
J j (B (div n) n* + grad*B (nn* - 1)) f ds + J J n* rot (Bn x C) 
and we obtain equation (A. 3) with div n = 



(Is 

J Jt. 

2 



We assume that E has no boundary and consequently, the term associated with 
r is null. 
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